Understanding the "anti-kick" in the merger of binary black holes 
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The generation of a large recoil velocity from the inspiral and merger of binary black holes represents one 
of the most exciting results of numerical-relativity calculations. While many aspects of this process have been 
investigated and explained, the "antikick", namely the sudden deceleration after the merger, has not yet found 
a simple explanation. We show that the antikick can be understood in terms of the radiation from a deformed 
black hole where the anisotropic curvature distribution on the horizon correlates with the direction and intensity 
of the recoil. Our analysis is focussed on Robinson-Trautman spacetimes and allows us to measure both the 
energies and momenta radiated in a gauge-invariant manner. At the same time, this simpler setup provides the 
qualitative and quantitative features of merging black holes, opening the way to a deeper understanding of the 
nonlinear dynamics of black-hole spacetimes. 
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Introduction. The merger of two black holes (BHs) is one of 
the most important sources of gravitational waves (GW) and 
it is generally accompanied by the recoil of the final BH as 
a result of anisotropic GW emission. While this scenario has 
been investigated for decades and first estimates have been 
made using approximated and semianalytical methods such as 
a particle approximation J2], post-Newtonian methods and 
the close-limit approximation (CLA) 0], it is only recently 
that accurate values for the recoil have been computed I5 H12I1 . 

Besides being a genuine nonlinear effect of general relativ- 
ity, the generation of a large recoil velocity during the merger 
of two BHs has a direct impact in astrophysics. Depending on 
its size and its variation with the mass ratio and spin, in fact, it 
can play an important role in the growth of supermassive BHs 
via mergers of galaxies and on the number of galaxies con- 
taining BHs 1130 . Numerical -relativity simulations of BHs in- 
spiralling on quasicircular orbits have already revealed many 
of the most important features of this process showing, for 
instance, that asymmetries in the mass can lead to recoil ve- 
locities Vk < 175 km/s lyllaL while asymmetries in the spins 
can lead respectively to < 450 km/s or v-^ < 4000 km/ s if 
the spins are aligned |8, 9, [Um or perpendicular to the orbital 
angular momentum QH-Hl (see lllql for a review). 

At the same time, however, there are a number of aspects of 
the nonlinear processes leading to the recoil that are far from 
being clarified even though interesting work has been recently 
carried out to investigate such aspects lfl7l - [l9Tl . One of these 
features, and possibly the most puzzling one, is the generic 
presence of an "antikick", namely, of one (or more) deceler- 
ations experienced by the recoiling BH. Such antikicks take 
place after a single apparent horizon (AH) has been found and 
have been reported in essentially all of the mergers simulated 
so far. 

This Letter is dedicated to elucidate the stages during which 
the antikick is generated and to provide a simple and qual- 
itative interpretation of the physics underlying this process. 
Our focus will be on the head-on collision of two nonspinning 
BHs with different mass and although this is the simplest sce- 
nario for a BH-merger, it contains many of the aspects that can 



be encountered in more generic conditions. Our qualitative 
picture will then be made quantitative and gauge-invariant by 
studying the logical equivalent of this process in the evolution 
of a Robinson-Trautman (RT) spacetime, with measurements 
of the recoil made at future null infinity. As commented be- 
low, the insight gained with RT spacetimes will be valuable to 
explain the antikick under generic conditions. 

The basic picture. Before discussing how to use the RT space- 
time to compute the antikick, it is useful to illustrate the basic 
BH physics leading to such process and for this we consider 
the collinear merger of two Schwarzschild BHs with unequal 
masses. This is shown in a schematic cartoon in Fig.[TJ where 
we have considered a reference frame in the centre of mass 
of the system and where the smaller black hole is initially 
on the positive z axis, while the larger one is on the negative 
axis. As the two BHs free-fall towards each other, the smaller 
one will move faster and will be more efficient in "forward- 
beaming" its GW emission 01. As a result, the linear mo- 
mentum will be radiated mostly downwards, thus leading to 
an upwards recoil of the BH binary [cf. stage (1) in Fig.[T|. At 
the merger the BH velocities will be larger and so will also be 
the anisotropic GW emission and the corresponding recoil of 
the system. However, when a single AH is formed comprising 
the two BHs, the curvature distribution on this 2-surface will 
be highly anisotropic, being higher in the upper hemisphere 
[cf. shading in stage (2) of Fig. [2. Because the newly formed 
BH will want to radiate all of its deviations away from the fi- 
nal Schwarzschild configuration, it will do so more effectively 
there where the curvature is larger, thus with a stronger emis- 
sion of GWs from the northern hemisphere. As a result, after 
the merger the linear momentum will be emitted mostly up- 
wards and this sudden sign change will lead to the antikick. 
The anisotropic GW emission will decay exponentially as the 
curvature gradients are erased and the BH will have reached 
its final and decelerated recoil velocity [cf. stage (3)]. 

Although this picture refers to a head-on collision, it is sup- 
ported by the findings in the CLA (where the direction of the 
ringdown kick is approximately opposite to that of the accu- 
mulated inspiral plus plunge kick) lfl7ll and it can be gener- 
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FIG. 1 : Cartoon of the generation of the antikick in the head-on col- 
lision of two unequal-mass Schwarzschild BHs. Initially the smaller 
BH moves faster and linear momentum is radiated mostly down- 
wards, thus leading to an upwards recoil of the system [stage (1)]. 
At the merger the curvature is higher in the upper hemisphere of the 
distorted BH (cf. shading) and linear momentum is radiated mostly 
upwards leading to the antikick [stage (2)]. The BH decelerates till a 
uniform curvature is restored on the horizon [stage (3)]. 

alized to a situation in which the BHs have different masses, 
different spins and are merging through an inspiral. Also in 
a more generic case, in fact, the newly formed AH will have 
a complicated but globally anisotropic distribution of the cur- 
vature, determining the direction (which is in general varying 
in time) along which the GWs will be emitted. Therefore we 
argue that the geometric properties in a dynamical horizon (of 
a black or white hole) determine its global dynamics. We next 
use the RT spacetime to validate this picture. 

The Robinson-Trautman spacetime. It is a class of vacuum 
solutions admitting a congruence of null geodesies which are 
twist and shear- free 02011 . with a future stationary horizon and 
a dynamical past (outer trapping) horizon l2lll (past AH here- 
after). A RT spacetime can thus be regarded as an isolated 
nonspherical white hole emitting GWs, where the evolution 
of the AH curvature-anisotropies and of the spacetime mo- 
mentum can be related unambiguously. The metric is 12211 
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where Q = Q(u, 0), u is the standard null coordinate, r 
is the affine parameter of the outgoing null geodesies, and 
fi = {9, 4>} are the angular coordinates on the unit sphere 
S 2 . Here M x is a constant and is related to the asymptotic 
mass, while the function K(u, fi) is the Gaussian curvature 
of the surface corresponding to r = 1 and u = constant, 
K(u, Si) = Q 2 (l + VqIiiQ), where is the Laplacian on 
S 2 . The Einstein equations then lead to 



d u Q(u,n) = -Q 3 ^lK{u,n)/(l2M a 



(2) 



Any regular initial data Q = Q(0, 0) will smoothly evolve 
according to (0 until it achieves a stationary configuration 
corresponding to a Schwarzschild BH at rest or moving with 
a constant speed l23tl . Equation (O implies the existence of 



the constant of motion A = J s2 dfl/Q 2 , which clearly repre- 
sents the area of the surface u,r = constant and can be used to 
normalise Q so that A = 47r. All the physically relevant infor- 
mation is contained in the function Q(u, ft), and this includes 
the gravitational radiation, which can be extracted by relating 
Q(u, SI) to the radiative part of the Riemann tensor B241 l25ll . 

The past AH radius R(u, ft) is given by the vanishing ex- 
pansion of the future ingoing null geodesies ll2"lll 
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The mass and momentum of the BH are computed at future 
null infinity using the Bondi 4-momentum lf22ll 
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withjr/"} = {l,sm#cos<^, sinflsin^, cos#}. Given smooth 
initial data, the spacetime will evolve to a stationary nonradia- 
tive solution which, in axisymmetry, has the form Q(oo, 9) = 
(lTvx)/Vl-v 2 , with x = cos 6» JH. The Bondi 4- 
momentum associated to Q(oo, 9) is 

{PM} Q = (Moo/a/1 - w 2 ) {1, 0, 0, ±v} , (5) 

so that the parameter v in Q(oo, 9) can be interpreted as the 
velocity of the Schwarzschild BH in the z-direction. 

One of the difficulties with RT spacetimes is the definition 
of physically meaningful initial data. Although we are more 
interested in a proof of principle than in a realistic configura- 
tion, we have adopted the prescription in l25ll 



Q(O,0) = Q n 
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which was interpreted to represent the final stages {i.e., after a 
common AH is formed) of a head-on collision of two boosted 
BHs with opposite velocities w and mass ratio q rf25tl . In prac- 
tice, to reproduce the situation shown in Fig. Q] we have set 
w < and taken q E [0, 1], but a more general class of initial 
data can be easily constructed. Note that Qo is normalized 
so that to A = 4tt and that in general the deformed BH will 
not be initially at rest. As a result, given the initial velocity 

v = P 3 (0)/P°(0), we perform a boost P° = A a p (v )pP 
3 

so that P (0) = by construction. The numerical solution of 
Eq. (|2]i with initial data (0 is performed as discussed in l22tl . 

Discussion. Figure [2] reports the typical evolution of a RT 
spacetime with the lower panel showing the evolution of the 
curvature of the past AH K An = 2M oc /R 3 (x) at the north 
(x = 1) and south pole (x = —1), and with the upper panel 
showing the evolution of the recoil velocity. Note that the 
two local curvatures are different initially, with the one in the 
upper hemisphere being larger than the one in the lower hemi- 
sphere (cf. Figure [TJ. However, as the gravitational radiation 
is emitted, this difference is erased. When this happens, the 
deceleration stops and the BH attains its asymptotic recoil ve- 
locity. The inset reports the curvature difference relative to the 
asymptotic Schwarzschild one, K Ali — 1, whose exponentially 
decaying behaviour is the one expected in a ringing BH. 
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FIG. 2: Typical evolution of a RT spacetime. Shown in the lower 
panel is the evolution of the curvature K AH at the north (x — 1) and 
south pole (x = —1). Shown in the upper panel is the evolution of 
the recoil, which stops decreasing when the curvature difference is 
erased by the radiation (dotted line). 

As mentioned before, that shown in Fig.|2]is a typical evolu- 
tion of a RT spacetime and is not specific of the initial data ((6]). 
By varying the values of w, in fact, it is possible to increase 
or decrease the final recoil, while a sign change in w simply 
inverts the curvature at the poles so that, for instance, initial 
data with w > would yield a BH accelerating in the posi- 
tive z-direction. Interestingly, it is even possible to fine-tune 
the parameter w so that the recoil produced for a RT space- 
time mimics the antikick produced by the quasicircular inspi- 
ral of nonspinning binaries. This is shown in Fig. [3j which 
reports the recoil as a function of the symmetric mass ratio 
v = q/(l + q) 2 , and where the dashed line refers to the anti- 
kick for the inspiral of nonspinning binaries in the CLA lfl7tl 
(the parameters chosen, i.e., w = —0.425 and m — 2 M, are 
those minimizing the differences). Considering that the two 
curves are related only logically and that the CLA one con- 
tains all the information about inspiralling BHs, including the 
orbital rotation, the match is surprisingly good. 

It is also suggestive to think that the curve in Fig.[3]is actu- 
ally composed of two different branches, one of which is char- 
acterized by large curvature gradients across the AH but small 
values of the curvature (this is the low-^ branch and is indi- 
cated with squares), while the other is characterized by small 
curvature gradients and large values of the curvature (this is 
the high-i^ branch and is indicated with circles). The same re- 
coil velocity can then be produced by two different values of 
v, for which the effects of large curvature gradients and small 
curvatures are the same as those produced by small curvature 
gradients but large curvatures. 

To go from this intuition to a mathematically well-defined 
measure we have computed the mass multipoles of the intrin- 
sic curvature of the initial data using the formalism developed 
in l26ll for dynamical horizons. Namely, we have calculated 



FIG. 3: Recoil velocity shown as a function of the symmetric mass 
ratio v when w = —0.425, with the dashed line refers to the antikick 
from the inspiral of nonspinning binaries in the CLA Jl7ll . Note 
that the curve can be thought of as being composed of two different 
branches. 

the mass moments as 

where P n (x) is the Legendre polynomial in terms of the coor- 
dinate x(0) which obeys dgx = - sin 9 R{8) 2 /(R 2 ah Q(9) 2 ), 
with R AH = y/A AH /(4:Tr) and x(0) = 1. Using these 
multipoles it is possible to construct an effective-curvature 
parameter K c g that represents a measure of the global cur- 
vature properties of the initial data and from which the re- 
coil depends in an injective way. Because this effective- 
curvature parameter has to contain the contribution from the 
even and odd multipoles, we have found that the expression 
A'off = M2I ^n=i ^2n+i/3 ra_1 |, reproduces exactly what is 
expected (note M% = to machine precision). 

This is shown in Fig. [4] which reports the recoil velocity as 
a function of K e g, As predicted, and in contrast with Fig. [3] 
the relation between the curvature and the recoil is now in- 
jective, with the maximum recoil velocity being given by the 
maximum value of K e g (see inset), and with the two branches 
coinciding. We do not expect the expression found here for 
K e g to be unique and indeed a more generic one will have 
to include also the mass-current multipoles to account for the 
spin contributions. However, lacking a rigorous mathematical 
guidance, our phenomenological K e g is a reasonable, intu- 
itive approximation. 

Conclusions. We have outlined a simple picture to explain 
the deceleration observed during the merger of binary BHs in 
terms of the dissipation of an anisotropic distribution of cur- 
vature on the horizon of the newly formed BH. We have ana- 
lyzed this picture for the head-on collision of two nonspinning 
BHs with unequal mass but its extension to generic systems is 
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FIG. 4: Recoil velocity shown as a function of the effective curvature. 
In contrast with Fig.[3] which uses the same symbols employed here, 
the relation between the curvature and the recoil is now injective. 

direct as the same features will be present also when includ- 



ing the spin and the orbital contributions: mass-current multi- 
poles will add (substract) in prograde (retrograde) orbits. The 
qualitative arguments made on the head-on collision have then 
been made quantitative by analyzing the gauge-independent 
dynamics of RT spacetimes. More specifically we have shown 
that the deceleration is associated to the radiation of curvature 
differences and persists as long as the gradients are not erased. 
Furthermore, the directionality of the recoil is dictated by the 
north-south curvature gradients and a one-to-one mapping be- 
tween the recoil and an effective curvature is possible. These 
results presented here can help in understanding some nonlin- 
ear aspects of curved spacetimes. 

Finally, an alternative interpretation of the recoil phe- 
nomenology can be given via the Landau-Lifshitz pseudoten- 
sor, where the recoil is given by the cancellation of large and 
opposite fluxes of momentum, part of which are "swallowed" 
by the BH lT27tl . While this is an interesting route, it relies 
on gauge-dependent measurements which may themselves be 
counterintuitive. 
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